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Finite temperature spatial qq correlation modes in the it and p channels are studied with 
the rainbow-ladder truncated quark Dyson-Schwinger equation and Bethe-Salpeter equa- 
tion in the Matsubara formalism. To retain the finite range of the effective interaction 
while facilitating summation over fermion Matsubara modes necessary to ensure conti- 
nuity at T = 0, a separable kernel is used. The model is fixed by T = properties and 
it implements dynamical chiral symmetry breaking and quark confinement. Above and 
below the deconfinement and chiral restoration transitions, we study M^-(T), f-w{T) and 
the 3-space transverse and longitudinal masses MT (T) and Mr (T) . For the p mode we 
also study the strong and electromagnetic partial widths to determine the T-dependence 
that is generated by the quark content. The p — * 7T7t width decreases sharply just above 
the chiral restoration transition leaving the p with only its narrow e + e~ width. We 
discuss improvements needed by this model, especially in regard to the high T behavior 
of the masses where we compare to lattice QCD simulations. 

PACS numbers: ll.10.St, ll.10.Wx, 12.38.Mh, 14.40.-n, 24.85. +p 

1. Introduction 

The connection between hot dense hadronic matter and a plasma of quarks and 
gluons is receiving increased attention with the advent of the relativistic heavy- 
ion collider fRHIC) at Brookhaven to complement previous investigations at the 
CERN SPS B. The plasma is expected to reveal itself through modified properties 
of hadronic reactions and their products. The di-lepton spectra has been of consid- 
erable interest □ as a relatively clean signal of how vector meson correlations and 
their decay channels and widths might be influenced by a hot and dense environ- 
ment. The question of how a vector meson strong decay, such as p — > mr, might 
respond as the temperature T or chemical potential ft crosses the critical phase 
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boundary for chiral restoration or quark deconfincmcnt requires that this process 
be studied at the quark-gluon level. To this end it is desirable to be able to describe 
quark deconfincmcnt and chiral restoration at finite T and p in a manner that can 
be extended to a variety of hadronic observables in a rapid and transparent way. 
The present model provides a simple framework for such investigations. 

At T = p = significant progress has been made within a continuum approach 
to modeling non-jperturbative QCD based on the truncated Dyson-Schwinger equa- 
tions (DSEs) Bala. An attractive feature of this approach is that dynamical chiral 
symmetry breaking and quark confinement can be embodied in the infrared struc- 
ture of the dressed gluon 2-point function which is constrained by a few chiral 
meson observables. Recent works have employed the ladder-rainbow truncation of 
the coupled quark DSE and qq Bethe-Salpeter equation (BSE) to produce success- 
ful descriptions of masses and decay constants of the light pseudoscalars ir and 
K U and the vectors p, (j>, and K* 13. These works have incorporated the one-loop 
renormalization group evolution of scale characteristic of QCD. With a few excep- 
tions U, applications to other hadron observables, including electromagnetic form 
factors a n d coup ling constants such as g p7 nr , have usually required the use of simpler 
models iffle. 

The finite T, p extension of realistic DSE/BSE models receives extra complica- 
tions due to the breaking of 0(4) symmetry and the dynamical coupling between 
Matsubara modes aO'La. These difficulties are compounded in the subsequent gen= 
eration of hadronic observables via straightforward adaptation of the approach □ 
found to be successful at T = p = 0. Studies of hadronic observables at finite T, p 
in DSE/BSE models have been restricted to simplifications such as extensions of 
the infrared-dominant (ID) model 113 in which the effective gluon propagator is re- 
stricted to an integrable singularity at zero momentum. Nevert heless s uch a model 
has proved capable of yielding qualitatively useful information tSlljIla. 

In this work we explore a separable Ansatz that can implement the essential 
qualitative features of DSE/BSE models at finite temperature. Separable represen- 
tations have previously been found capable of an efficient modeling. of the effective 
qq interaction in the infrared domain for T = meson observables liScj. A previous 
implementation at T > employ 
a quark confinement mechanism I 



an instantaneous separable interaction without 
Here the approach is covariant and the simple 
separable interaction absolutely confines quarks at T = 0. The few parameters are 
fixed by_7r and p/ui properties. The approach is a simplification of one developed 
earlier 123 that was found to be quite successful for the light meson spectrum at 
T = 0. The confining mechanism is an infrared enhancement in the quark-quark 
interaction that is strong enough to remove the possibility of a mass shell pole in 
the quark propagator for real p 2 . In the simple implementation used here, it is 
particularly transparent that sufficiently high temperature will necessarily restore a 
quark mass-shell pole and there will be a deconfinement transition. This model also 
implements low temperature dynamical chiral symmetry breaking and it preserves 
the Goldstone theorem in that the generated ir is massless in the chiral limit. The 
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solutions of the BSE for the 7r and p modes are particularly simple and are used to 
study the T-dependence of the meson masses and decay constants in the presence 
of both deconfincment and chiral restoration mechanisms. Some preliminary results 
have been previously discussed for T > E3 and for T, p > Ej. 

In Sec. II the T = separable model is introduced, the DSE solutions for the 
dressed quark propagator are developed and the quark confinement property is de- 
scribed. The BSE solutions for 7r and p are also treated there and the corresponding 
decay constants are defined. The extension to T > is considered in Sec. Ill and the 
chiral symmetry restoration and deconfinement phenomena are discussed. The spa- 
tial 7r and p modes are considered there. The Gell-Mann-Oakes-Renner (GMOR) 
relation and its generalization are used to evaluate the performance of this model 
in respecting chiral symmetry constraints. The widths of the transverse p from 
p — > e + e~ and p —f ittt are also considered in Sec. III. The high T behavior of the 
model is considered in Sec. IV by comparison of masses with results from lattice sim- 
ulations. A discussion follows in Sec. V. An Appendix details the high T behavior 
of the ID model as a reference. 

2. Confining separable Dyson-Schwinger equation model 

Mesons can be described as qq bound states using the Bethe-Salpeter equation. 
In the ladder truncation, this equation readsffl 

- X(P 2 )T(p,P) = ±J j£jLDf v (p- qh ll S(q + nq,P)S(q-) lv , (1) 

where P is the total momentum, q± = q ± P/2, and D^l(k) an "effective gluon 
propagator". The meson mass is identified from X(P 2 = — M 2 ) = 1. In conjunction 
with the rainbow truncation for the quark DSE 

Sip)' 1 = Z 2 z 1 -p + Z 2 m Q + ^J^g 2 Df u (p~ q ) 1 ,S( q ) lu . (2) 

this equation forms the basis for the DSE approach to meson physics US. 

Recent studies have employed D c ^(k) = G{k 2 )D^° c (k) where D f ™ e {k) is the free 
gluon propagator in Landau gauge. The effective coupling Q (k 2 ) is given by the one- 
loop perturbative form of the QCD running coupling in the ultraviolet while the 
phenomenological infrared form is chosen to reproduce the pion and kaon masses 
and decay constants eL This model can succesfully describe the vector meson masses 
and dominant decays 0, as well as the pion and kaon electromagnetic form factors B, 
without new parameters. 

The direct extension of such an approach to accommodate finite temperature E3 
and baryon density entails a significant increase in complexity. In the Matsubara 
formalism, the number of coupled equations represented by Eqs. ([!]) and (^) scales 
up with the number of fermion Matsubara modes included. For studies near and 

t\Ve use a Euclidean space formulation, with {7^,7^} = 2<5 MJ/ , 7J = 7^ and a ■ b = Oj&j. 
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above the transition, T > 100 MeV, about 10 such modes appear adequate El The 
appropriate number can be ~ 10 3 if reasonable continuity with T = results is to 
be verified. Meson qq modes at T > have often been studied within the Nambu- 
Jona-Lasinio model where the contact nature of the effective interaction allows 
decoupling of theMatsubara modes and also analytic methods of summation, see 
for example, Ref. E9 and references therein. However the lack of quark confinement 
in that model leads to unphysical thresholds for qq dissociation. 

We consider here a simple separable interaction that has a finite range, ac- 
commodates quark confinement, and facilitates a decoupling offermion Matsubara 
modes. We base our approach on a confining separable model E3 previously found 
to be successful at T = and defined by D c ^,(p — q) — > 5^ D(p 2 , q 2 ,p ■ q) with 

D(p 2 ,q 2 ,p-q) = D f (p 2 )fo(q 2 ) + D 1 h{p 2 )(p ■ q)h{q 2 ) . (3) 

Here a Feynman-like gauge is chosen for phenomenological simplicity. This is a 
rank- 2 interaction with two strength parameters Dq,Di, and corresponding form 
factors fi{p 2 )- The choice for these quantities is constrained by consideration of the 
resulting solution of the DSE for the quark propagator in the rainbow approxima-_ 
tion. For the amplitudes defined by S{p) — [i^A{p 2 ) + B(p 2 ) + mo] -1 this produces! 

3 7 (2tt) 4 q 2 A 2 (q 2 ) + [B{q 2 ) + m ] 

LV ' J 3 7 (2tt) 4 KP ' p q ' q ^A 2 (q 2 ) + [B(q 2 )+m a } 2 W 

We note that if terms of higher order in p ■ q were to be included in ^n. (^), 
they would make no contribution to Eqs. (||) and ^ for the DSE solution E3. The 
solution for B(p 2 ) is determined only by the Do term, and the solution for A{p 2 ) — 1 
is determined only by the D\ term. Eq. (|^) produces solutions having the form 

B(p 2 ) = b f (p 2 ) , A(p 2 ) = 1 + a h{p 2 ) , (6) 

and Eqs. (Q) and (JsJ) reduce to nonlinear equations for the constants b and a. 

2.1. Confinement and Dynamical Chiral Symmetry Breaking 

If there are no poles in the quark propagator S(p) for real timelike p 2 then there 
is no physical quark mass shell, quarks cannot propagate free of interactions, and the 
description of hadronic processes will not be hindered by spurious quark production 
thresholds. This is sufficient but not necessary for quark confinement and it remains 
a viable possibility for how quark confinement is realized 00. The more general 
phenomena of confinement of colored multi-quark states and the non-existence of 
S-matrix elements connecting them to hadronic states are more subtle topics that 



§ We choose the interaction form factors such that they provide sufficient ultraviolet suppression. 
Therefore no renormalization is needed and Zi = 1 . 
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do not concern us hereu. A nontrivial solution for B(p 2 ) in the chiral limit (mo = 0) 
signals dynamical chiral symmetry breaking. There is a connection between quark 
confinement as the lack of a quark mass shell and the existence of a strong quark 
mass function in the infrared through dynamical chiral symmetry breaking. This 
connection has proved to be empirically successful for the description of ground 
state light quark mesons and their form factors and decays BEil. In the present 
separable model the strength b = B(0), which is generated by solution of Eqs. (Q) 
and (^), controls both confinement and dynamical chiral symmetry breaking. 

The propagator is confining if m 2 (p 2 ) ^ —p 2 for real p 2 where the quark mass 
function is m{p 2 ) — (B(p 2 ) + uiq) / A(p 2 ) . With an exponential form factor fo(p 2 ) — 
exp (— p 2 /Aq), this condition is most transparent in the case of a rank-1 separable 
model where D\ = and A(p 2 ) — 1, i.e., a = 0. In the chiral limit, the model is 
confining if Do is strong enough to make &/Ao > 1/ V2e and this finite b also signals 
dynamical chiral symmetry breaking. Using Ao ~ 0.6 — 0.8 GeV as a typical range 
for a quark mass function m{p 2 ), both confinement and dynamical chiral symmetry 
breaking will be compatible with m(p = 0) > 0.3 GeV, an empirically viable value. 
The above qualitative properties will also hold in the case of a rank-2 model if the 
form factor ranges satisfy Ai > Ao and this is compatible with empirical findings 
that the amplitude A(p 2 ) typically has a larger momentum range than B(p 2 ). At 
finite temperature, the strength b(T) for the quark mass function will decrease with 
T so that this model can be expected to have a deconfincmcnt transition at or before 
the chiral restoration transition associated with b(T) — > 0. 

It is found that the simple choice fi(p 2 ) — exp (— p 2 /Af) produces numerical 
solutions that describe infrared properties of ir and u> mesons very well and generate 
an empirically acceptable chiral condensate. At the same time the produced quark 
propagator is found to be confining and the infrared strength and shape of the quark 
amplitudes Alp 2 ) and B(p 2 ) are in qualitative agreement with the results of typical 
DSE studies El. We use the exponential form factors as a minimal way to preserve 
these properties while realizing that the ultraviolet suppression is much greater 
than the power law fall-off (with logarithmic corrections) known from asymptotic 
QCD. Most of our investigation centers on physics below and in the vicinity of the 
transition. We use the high T behavior of masses in comparision with lattice results 
to discuss improvements appropriate for future work. 

2.2. 7r and p bound states 

With the separable interaction of Eq. (||) , thii allowed form of the solution of 
Eq. (0) for the pion BS amplitude TT„{q;P) is E9 

I\,(g;P) = 75 (iE„{P 2 ) + PF^P 2 )) fo(q 2 ) , (7) 

which contains the two dominant covariants from the set of four general covariants. 
The q dependence is described only by the first form factor fo(q 2 )- The second 
term /i of the interaction can contribute only indirectly via the quark propagators. 
The 7T BSE, Eq. (@), 

becomes a 2 x 2 matrix eigenvalue problem JC(P 2 )j = X(P 2 )f 
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where the eigenvector is / = (E^ , ) . The kernel is 



K ii {P i ) = - 



4Dn 



tr* 



(8) 



where the it covariants are t = (275,75 P) and we have also introduced t = (175, 
—75 P/2P 2 ). We note that the separable model produces the same q 2 shape for 
both amplitudes F^ and E v ; the shape is that of the quark amplitude B(q 2 ). For the 
general amplitude E^{q; P) this is the correct shape fn the chiral limit; for physical 
quark masses it is still a very good approximation In general, the amplitude 
F v {q\P) does not have that shape; it is in fact linked with A(q 2 ) — 1 through the 
axial vector Ward-Takahashi identity (AV-WTI) c£l. However Goldstone's theorem, 
the key consequence of the AV-WTI, is preserved by the present separable model; 
in the chiral limit, whenever a nontrivial DSE solution for B(p 2 ) exists, there will 
be a massless 7r solution to Eq. 

A simple illustrative truncation is obtained by setting F^iP 2 ) = for then 
Eq. (0) reduces to an expression for the eigenvalue which is 



K(P 2 ) = 



I6D0 



(27T) 



f/oV 



P< 



>v u v 



'S U S 



(9) 



where the quark propagator amplitudes employed here are defined by S(p) = — i$ 
°v(p 2 )+ o"s(p 2 )- Since A = 1 in rank-1, we omit the amplitude F^ in that case. 
Thus it is Eq. (|9|) that we use for the 7r calculation with a rank-1 interaction. 

The vector mesons p and lo are degenerate in the ladder approximation. We 
shall deal with the isovector o. For the BS amplitude fr^(g; P) there are in general 
eight transverse covariants and the dominant one is yE{P) = T fJ , u (P) r y u where 



J fJ,V 



PpPv/P ■ The solution of the rank-1 model contains only that 



term, that is, T^(q; P) 

by 



7u (P)fo(Q )Fp{P )■ The corresponding eigenvalue is given 



A P (P 2 



8D f d 4 q 



p2 

4 



2q 2 



(1 



Vv a v 



CT S a S 



(10) 



where z = q ■ P. For the rank-2 separable interaction, there are two other covariants 
besides 7J that will appear in the solution of the BSE El. However it has been found 
in such a model that these subleading vector covariants make only a few percent 
contribution to the vector mass E3 and the associated g p7r7r 0- In this present work 
we will ignore the subdominant covariants and employ 7J as the only covariant for 
the vector meson. The differences in the vector mass obtained from Eq. (^) in 
rank-1 and rank-2 will be due to differences in the quark propagator in each case. 
The normalization condition for the ir BS amplitude can be expressed as 



2P„ 



d 
dP L 



2N n tr. 



d 4 q 
(2nY 



t 7r (q;-K)S(q + )T 7r (q;K)S(q. 



(11) 



P 2 =K 2 = 



-Ml 
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Here f(q;K) is the charge conjugate amplitude [C~ 1 T(— q, K)C\ l where C =7274 
and t denotes matrix transpose. This defines a normalization constant N n via 
E 7r (P 2 = —M 2 ) fo(q 2 ) = B(q 2 )/N.K- The pion decay constant f n can be expressed 
as the loop integral 

UP„Sii = (0|g^7„7Mki(i , )> 

= 8 l0 N c tr s J — 757m r„(g; P) S(g_) . (12) 

Note that Eq. ([l2|) is the exact expression for /„. except for the absence of the renor- 
malization constant Z2 which is unity in the present model. The decay constants 
thus test the quality of the infrared behavior of the quark propagator and the BS 
amplitudes. Similar expressions exist for the normalization of the vector meson BS 
amplitude and the coupling between a vector meson and a photon 0O. 

Table 1. Calculated properties at T = for the u/d quark mesons it and p along with related quark 
properties and parameters for the rank-1 and rank-2 separable models. The quoted experimental 
values for both the quark condensate and the current u/d quark mass mo are appropriate to a 
renormalization scale of fi ~ 1 GeV. The scale for the model calculations may be considered to be 
set by Aq. 





Experiment 


rank-1 


rank-2 


- <99> u 


(0.236 GeV) J 


0.248 


0.203 


mo 


5 - 10 MeV 


6.6 


5.3 


m(p 2 = 0) 


~0.350 GeV 


0.685 


0.405 


Mir 


0.1385 GeV 


0.140 


0.139 




0.093 GeV 


0.104 


0.093 




1.0 


0.987 


0.740 




0.770/0.783 GeV 


0.783 


0.784 


9p 


5.04 


5.04 


6.38 


r p o_ e + e - 


6.77 kcV 


6.76 


4.22 




6.05 


5.71 


7.22 


— >7T7r 


151 MeV 


137 


221 


Parameters 




A)A.o 




128.0 


260.0 


Ao 




0.687 GeV 


0.638 GeV 


DiAf 







130.0 


A1/A0 






1.21 



In Table [I] the results for the ground state 7r and p/uu mesons as well as related 
quantities are shown for both rank-1 and rank-2 versions of the model along with 
the values of the employed parameters. We consider the experimental u> mass to 
be the appropriate value for comparison with the vector result in ladder approx- 
imation in order to allow for subsequent preferential lowering of M p due to pion 
loop dressing cQ. The range parameter Ao is used to set the mass scale to repro- 
duce the experimental M w for rank-1 and f w for rank-2. The obtained dynamical 
quark mass function at p 2 = is m(p 2 = 0) = 0.685 GeV for rank-1. For rank-2 
we obtain A(p 2 = 0) = 1.94 and m(p 2 = 0) = 0.405 GeV. In both cases the dressed 
quark propagator is confining. The strengths obtained for the mass function are 
consistent with results from recent DSE solutions &Q. Also shown in Table |] are 
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the results obtained for the electromagnetic (g p ) and strong decays {g P Tr-n) of the p, 
which compare reasonably well with experiments. 

The difference between N v and f„ shown in Table |l| simply reflects the fact that 
the AV-WTI cannot be exactly satisfied within a separable model. (This is also 
the case in the Nambu-Jona-Lasinio model due to the required cut-off and in any 
approach that does not use a translationally invariant interaction or_a sufficiently 
complete set of Dirac covariants for and the vertex amplitudes cB.) A trans- 
lationally invariant model that makes the simplifying assumption A(p 2 ) = 1 can 
achieve a fortuitious agreement u and this dominates the rank-1 result for N n /f„. 
Rather than use f n to define N w , we calculate N„ from its definition in Eq. ( pT| ) so 
that, in subsequent studies of processes such as p — > 7T7t, the pion state is physically 
and consistently normed within the model. 

3. Finite Temperature Extension 

The extension of the separable model studies to T ^ is systematically accom- 
plished by transcription of the Euclidean quark 4-momentum via (u> n , q), 
where uj n = (2n + 1)tcT are the discrete Matsubara frequencies. The effective qq 
interaction will automatically decrease with increasing T without the introduction 
of an explicit T-dependence which would require new parameters. We investigate 
the resulting behavior of the 7r and p meson modes and decays in the presence of 
deconfincmcnt and chiral restoration. 

3.1. Chiral symmetry restoration and deconfinement 

The result of the DSE solution for the dressed quark propagator now becomes 

S~ x {p n , T) = ij -pA(pl T) + i l4 u n C(p 2 n , T) + B(p 2 n , T)+m , (13) 

where p 2 n = Lo 2 n + p 2 and there are now three amplitudes due to the loss of 0(4) 
symmetry. The solutions have the form B — b(T) fo(p 2 l ) , A = 1 + a(T) fi(p 2 l ) , and 
C = 1 + c(T) fi(p 2 l ) and the DSE becomes a set of three non-linear equations for 
b(T), a(T) and c(T). The explicit form is 

a(T) = ^Tj2 J fxipDfil + a^Mpl^d-'iplT) , (14) 
c(T) = ^ T E / ^fM)"i[l + c{T)h(pl)] d-\p 2 n ,T) , (15) 

b{T)= lGDo T J2 [ ^_ /o (p2 )K + 6(T)/o(p 2 )]rf -l (p 2 r); (16) 

where d(p 2 nl T) is given by 



d(p 2 n , T) = p 2 A 2 (p 2 n , T) + u 2 n C 2 (p 2 n , T) + [m + B(p 2 n , T)f 



(17) 
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T [GeV] 



Figure 1: T-dependence of quark self-energy amplitudes at p = from solution of 
the DSE. 



It is also useful to introduce the equivalent representation 

S(p n ,T) = -tf ■ p a A (pl,T) -ij 4 u n <Jc(pI,T) + <T B (pl,T) , (18) 

where a a — A/d, oc = C/d and <jb = (B + mo)/d. 

The T-dependence of the solutions for A, B and C at p 1 = for the rank-2 
model is displayed in Fig. [l]. The results shown are for the lowest Matsubara mode 
(n = 0) which provides the leading behavior as T is increased. The chiral restoration 
critical temperature T c is identified from the vanishing of the chiral limit amplitude 
B Q (p = 0,T) as shown. We find T c = 121 MeV for the rank-2 model. Below T c , A 
and C are relatively constant, 0(4) symmetry is approximately manifest, and the 
main effect is an almost constant quark wave function renormalization via 1/(7; the 
central feature shared by both rank-1 and rank-2 models is a rapidly decreasing 
mass function. Above T c , there remains a significant temperature range where the 
self-interaction effects are strong, both A and C are considerably enhanced above 
their perturbative values, and the breaking of 0(4) symmetry is manifest. The 
present model thus captures the qualitative T-denendence observed for the dressed 
quark propagator in studies of the quark DSELLjO The rank-1 limit has A = C = 1 
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for all T and the behavior of B(p = 0,T) is similar to that in Fig. [j] except T c = 
146 MeV is obtained. 

The order parameter for chiral restoration, the chiral quark condensate, can be 
obtained from the chiral limit quark propagator as (qq)° = —N c tx s So(x,x). Here 
this produces 

where bo and do are obtained from the chiral limit solution of the DSE. Both bo(T) 
and (qq}° vanish sharply as (1 — T/T C )P with the critical exponent having the mean 
field value (3 = 1/2 in agreement with other rainbow DSE studies Eil. 

The deconfinement temperature is found by a search for a propagator pole 
(a zero of the function c?(pg,T)) on the real p 2 axis. We find Td — T c = 146 MeV 
for the rank-1 model and Td = 0.9 T c = 105 MeV for the rank-2 model. Only about 
15% variation in these transition temperatures can be achieved by variation of 
the model parameters while retaining a reasonable description of the observables 
shown in Table [j]. Since all dynamical information concerning the rank-1 model 
is contained in one function B{p 2 ll T), one may expect to find Td = T C in that 
case. For rank-2, the dynamical information is contained in three functions and the 
result Td ^ T c is not surprising. It is however possible that the separable form of 
interaction used here might miss some dynamical correlations between A, B and C 
that would otherwise produce Td = T c . 

3.2. Spatial it correlations at T ^ 

At T = the mass-shell condition for a meson as a qq bound state of the BSE is 
equivalent to the appearance of a pole in the qq scattering amplitude as a function 
of P 2 . At T 7^ in the Matsubara formalism, the 0(4) symmetry is broken by 
the heat bath and we have P — > (f2 m , P) where f2 TO = Im-nT . Bound states and the 
poles they generate in propagators may be investigated through polarization tensors, 
correlators or Bethe-Salpeter eigenvalues. This pole structure is characterized by 
information at discrete points Q m on the imaginary energy axis and at a continuum 
of 3-momenta. Analytic continuation for construction of real-time Green's functions 
(and related propagation properties) has been well-studied 123. An unambiguous 
result is obtained by requiring that the continuationyield a function that is bounded 
at complex infinity and analytic off the real axis E3. One may search for poles as 
a function of P 2 thus identifying the so-called spatial or screening masses for each 
Matsubara mode. These serve as one particular characterization of the propagator 
and the T > bound states. 

In the present context the eigenvalues of the meson BSE become A(P 2 ) — > 
A(ri^, P 2 ;T). The temporal meson masses identified by zeros of 1 — A(£l 2 , 0; T) will 
be different in general from the spatial masses identified by zeros of 1 — A(0, P 2 ; T). 
They are however identical at T = and an approximate degeneracy can be ex- 
pected to extend over the finite T domain where the 0(4) symmetry is not strongly 
broken. From Fig. |l| one may reasonably expect this domain to extend up to 
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T ~ 80 — 100 MeV. At and above the transition, temporal and spatial masses can 
be expected to emphasize different aspects of the bound state modes. In this work 
we explore the T-dependence of the lowest spatial masses in the 7r and p channels. 

In the 7r channel, the correlator Ilfx) = (T J ps (x) J ps (0)} of two currents J ps (x) — 
q(x)"f§q(x), after transformation to momentum space and extension to T > via 
the Matsubara formalism, can be expressed as 

U(n 2 ,P 2 )=N c T^2tT s J ^ l5 S(q n+ )T ps ( qn - i n,P)S(q n -), (20) 

where q n ± — q n ± P/2 and T ps is the pseudoscalar vertex which satisfies the in- 
homogeneous version of the pion BSE. The quark propagators appearing here are 
dressed according to solution of the DSE. When the ladder approximation is em- 
ployed for T psi along with the present separable approximation to the BSE kernel, 
the correlation function can be expressed as 

n(n 2 ,p 2 ) = rt°\n 2 1 P 2 )~^N c L l (n 2 ,P 2 ) [i-jc(n a , J 32 )],-. 1 Un 2 ,P 2 ) . (21) 

3 1 j 

Here II^ ) is the contribution to Eq. (|2(]) arising from the zeroth order contribution 
(75) to the pseudoscalar vertex T ps . The second term of Eq. (^T|) sums the inter- 
action terms and factorizes due the separability of the effective interaction. The 
kernel JC involves the T > extension of the kernel of the tt separable BSE given 
previously in Eq. (|8|). The loop integrals for the numerator are given by 

L i (n 2 ,P 2 )=T j2^sJ -^f (q 2 n )^ S(q n+ )US(q n -), (22) 

and 

L J (n 2 ,P 2 )=T^sJ ^h{ql)hS{qn + )l,S{q^). (23) 

With an eigenvector representation of the BSE kernel /C, Eq. ( |2l| ) develops a 
denominator 1 — A Tr (il 2 , P 2 ; T). There is a pole in the correlator associated with 
the spatial mode solution to the homogeneous BSE identified from 

1 - \ w (0,P 2 ;T) = Z- 1 {P 2 ,T)[P 2 + M 2 (T)} - . (24) 

The masses so identified are spatial screening masses of the lowest mode associated 
with the 3-space asymptotic behavior H(x) ~ exp(— Mx). The identification of 
spatial masses by location of a pole is equivalent, by Fourier transformation, to the 
method of large spatial separation of sources used in lattice QCD. 

The general form of the finite T pion BS amplitude allowed by the separable 
model is 



T„(q n ; P m ) = 75 + 74 ^ m F n (P^) + 7 • PF r (P*j) h{q 2 n ) ■ (25) 
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The separable BSE becomes a 3 x 3 matrix eigenvalue problem with a kernel that is 
a generalization of Eq. (||). In the limit f2 m — > 0, as is required for the spatial mode 
of interest here, the amplitude F n — ^l m F v is trivially zero. The pseudovector 
amplitude F„ is significantly different from zero below T c , but decreases rapidly 
above the transition. In the chiral limit, it vanishes identically at and beyond the 
transition; above the transition, only is survives. 

The result for the 7r mass is displayed in Fig. |]for rank-2; the results for the rank- 
1 model are similar. In both models, M^{T) is seen to be only weakly T-dependent 
until near T c where a sharp rise begins. These qualitative features of the response of 
the pion mode with T agree with the results deduced from the DSE in Ref. 113 and 
also with the more detailed study of Ref. E3 using a ladder-rainbow truncation of the 
DSE/BSE system that preserves the one-loop renormalization group properties of 
QCD. Evidently the detailed character of the effective interaction in the perturbative 
region does not dominate at the level of the present qualitative investigation. 
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Figure 2: Chiral symmetry in the rank-2 model. T-dependence of M v and quantities 
involved in the exact pion mass relation and the GMOR relation. 



3.3. Chiral symmetry and it mass relation 

To explore the extent to which the model respects the detailed constraints from 
chiral symmetry, we investigate the exact QCD pseudoscalar mass relation cB which, 
after extension to the spatial mode at T > 0, is 



Ml(T)U{T) = 2m r P {T) 



(26) 
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Here rp, the residue at the pion pole in the pseudoscalar vertex, is given by the pseu- 
doscalar projection of the pion wavefunction onto zero quark-antiquark separation, 
that is 

ir P (T)=N c Tj2^ s f 44al5S(q n + ^)T^q n -P)S(q n -^). (27) 

The generalization of Eq. ( |i2|) for /„. to finite temperature in the case of the spatial 
pion mode is 

C d?a P P 

PiU(T) = N c T^tTsJ 757* S(q n + -) T„{q n ;P) S(q n - -) . (28) 



Eqs. ( J2 7| ) and (|2§|) are exact expressions for rp(T) and f^iT) except for the absence 
of the renormalization constants which are trivially equal to one in this separable 
model. The relation in Eq. (^6|) is a consequence of the pion pole structure of the 
isovector axial Ward identity whichlinks the quark propagator, the pseudoscalar 
vertex and the axial vector vertex cB. In the chiral limit, rp — > —(qq)°/f® and 
Eq.(^6|), for small mass, produces the Gell-Mann-Oakes-Renner (GMOR) relation. 
The exact mass relation, Eq. (^6|), can only be satisfied approximately when the 
various quantities are obtained in a manner that does not preserve axial Ward 
identity. The error can be used to assess the reliability of the present approach to 
modeling the behavior of the pion spatial mode as the temperature is varied. 

Our findings in the case of the rank-2 model are displayed in Fig. |[ There 
the solid line represents rp(T) calculated from the quark loop integral in Eq. (27); 
the dotted line represents rp constructed from the other quantities in Eq. (2£). 
The exact mass relation, Eq. (^6|) , is violated by about 25% almost independent of 
temperature, even above T c . The rank-1 model satisfies this relation to within 1%, 
due to the special but unrealistic case A = 1 = C . We have also investigated the 
(approximate) GMOR relation within the present model. The quantity —(qq)°/N®, 
displayed in Fig. ||, is the chiral limit of rp in this model. If all covariants for 
the pion were retained and the axial vector Ward identity were obeyed, one would 
have iV° = f® in the chiral limit cB. If the GMOR relation were exactly obeyed, 
the long-dashed line representing (qq)°/N® would coincide with the dotted line 
representing M% f^/2'mQ. These features are temperature-independent until about 
0.9 T c , consistent with an earlier study of low energy theorems at T > within a 
three-space, non-confining separable interaction model H. Close to T c , the GMOR 
relation breaks down, and above T c , it is no longer well-defined. 

It should be noted that /°, N® and (qq)° are equivalent order parameters for the 
critical behavior near T c and have weak T-dependence below T c . A consequence is 
that M^/ T , rp and (qq)°/N® are almost T-independent and so are the estimated 
errors for the two mass relations linking these quantities. We have employed the 
physical /„■ defined at non-zero current quark mass, and this does not vanish at T c 
but continuously decreases. 
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3.4. Spatial p correlations at T ^ 

The T — transverse vector meson, that we have described by the covariant 
7J, splits for T > into 3-space longitudinal and transverse modes. For the spatial 
modes characterized by P = (0, P) the BS amplitudes are 

T^ L \q n ;P) = ^4 74/ofe 2 l )F p(i) (P 2 ) , (29) 

and 

rf T) (q„l?) = (li - fo(ql)F p(T) (P 2 ) . (30) 

The T-dependence of the corresponding masses is displayed in Fig. ^ for the rank-2 
model. These modes are effectively degenerate and T-independent until about T c /2 
where the breaking of 0(4) invariance becomes significant. The qualitative features 
M^{T) > Mj(T) and Mj(T) sa const for T <T C seen here in the present context 
of a finite range interaction have previously been noted within the limiting case 
of the zero momentum range ID model 113. This latter model was not applied for 
T > T c . We discontinue the present study of the longitudinal mode at T ~ 180 
MeV where it becomes unstable to qq dissociation. The transverse mode continues 
to be below the spatial qq threshold for the temperature range displayed. 

The T — expression EB for the electromagnetic coupling constant g p has a 
straightforward extension to T > for a transverse spatial p° mode. Use of the 
dm = solution described above yields 



9 P {T) 



■ £ trs / 7i s{qn + ? } T * m {qn] p) s[qn ~ ? } ' (31) 



With summation over enough Matsubara modes for convergence, the produced 
g p {T) tends smoothly to the previously determined T — result. The result over a 
temperature range that extends just beyond T c is displayed in Fig. [|. There it is 
confirmed that there is very little T-dependence below about 0.9 T c where there is 
approximate 0(4) symmetry as evident in the quark propagator behavior in Fig. ^. 



The impulse approximation for the pirir vertex aliil, after extension to T > for 
spatial modes characterized by Q = (0, Q) for the p and P = (0, P) for the relative 
7T7T momentum, takes the form 

&„(P,Q) - Pu (T) 

= -2N c T^2tr s J ^T„(k n+ ;-P + )S(q n+ -)T^ L Hq n+ ;Q) 

xS(q n++ ) TAkn-;P-) S{q n J) . (32) 

Use of Eq. ( f^ ) immediately shows that the longitudinal p mode cannot couple to 
7T7T. The temperature dependence obtained for the transverse p coupling constant 
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Figure 3: T-dependence of the longitudinal p mass (dot-dashed line) and the trans- 
verse p mass (solid line). Also shown is 2M„ indicating that the strong decay 
channel becomes inaccessible within 20-30 MeV beyond T c . Results are from the 
rank- 2 model. 



9pirir(T) is displayed in Fig. ||. Again one observes continuity with the T = result 
and a very weak T-dependence until about 0.8T C . Around T c , the coupling constant 
decreases significantly. 

For both interactions of the p mode, and below T c , one expects qualitatively 
similar behavior from spatial modes associated with higher meson Matsubara fre- 
quencies fl m 7^ except that the effect of 0(4) symmetry breaking will be evident 
earlier. Above T c , each 2-7rT increment to the meson Matsubara frequency adds 
significantly to the quark effective mass in integrals like Eq. (|3l]). Since dynam- 
ical chiral symmetry breaking is now absent, the quark propagators are of Dirac 
vector character, the meson Matsubara frequencies are the largest mass scale in 
the system, and perturbative behavior becomes increasingly dominant. The various 
spatial modes from fl m ^ are characterized by masses much greater than that of 
the lowest mode considered here. We anticipate that this lowest mode characterizes 
the qualitative behavior of the physical decay processes p° — > e + e~ and p — ► hit. 
Certainly, a high temperature limit in which some modes vanish and others diverge 
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Figure 4: T-dependence of the p electromagnetic decay constant g p (dashed line) 
and the strong coupling constant g p7T T, (solid line). Results are from the rank-2 
model. 



would indicate a non-analytic behavior in the reconstructed physical amplitude that 
is physically untenable. We therefore estimate the relevant decay widths by com- 
bining the present results for the coupling constants with the relevant phase space 
factors generated also from the lowest spatial mass modes. 
This leads to the electromagnetic decay width 



r p o^ e + e _ (T) 



4WMj(T) 



3 9 P {Tf 



while the corresponding strong decay width is 



r(T) 



5p _(T) 2 Mj(T) 



An 



12 



4M 7r (r) 2 

Mj{Tf 



3/2 



(33) 



(34) 



The T-dependence estimated in this way for the decay widths is due to the re- 
sponse of the quark substructure to the heat bath, particularly the restoration of 
chiral symmetry. The results are displayed in Fig. ^. The contrast between the 
behavior of the electromagnetic and strong widths near and just above T c should 
be a more robust finding than the details of the individual processes. The strong 
width decreases rapidly and vanishes just above T c while the electromagnetic width 
remains within 20% of the T = value. Part of the strong decrease of the intrinsic 
7T7T width of the transverse p is due to the decrease in the coupling constant, however 
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Figure 5: T-dependence of the transverse p partial widths due to electromagnetic 
e + e~ decay (dashed line) and strong tttt decay (solid line) corresponding to the 
coupling constants shown in Fig. ^. 



the dominant effect is the T-dependence of the last factor in Eq. ( |34| ) . As displayed 
in Fig. |, 2M 7r (T) rises faster with T than does Mj(T) until at T = 1.17 T c we have 
Mj = 2M n . Beyond this point, the phase space factor vanishes and the strong de- 
cay p T — > 7T7T is blocked. This suggests that the total p T width of 151 MeV at T = 
decreases by about 50% near T — T c and drops sharply to the electromagnetic value 
of about 6 keV by T = 1.17 T c . 

This narrowing of the intrinsic decay width of the vector meson mode in the 
heat bath is a mechanism that is distinct from the collisional broadening effect 13 
from the many-hadron environment. The present work indicates that there is a 
non-trivial T-dependence to intrinsic coupling constants such as (? p7r7r and decay 
phase space. The intrinsic effect tends to significantly decrease the decay width; 
the many-hadron medium effects have the opposite influence. Phenomenological 
forms for the T-dependence of M p (T) and r p o_ W7r (T) have often been explored in 
studies of medium effects in heavy-ion collisions. For example, both an increase 
in the width of the form 1^/(1 — T 2 /T 2 ), and a p mass decreasing by 50% at T c 
ve been explored in an effort to understand the heavy-ion dilepton spectrum 
A coordinated approach is called for in which hadronic collisional broadening 
mechanisms are built upon intrinsic coupling constants that respect the temperature 
and density dependence of the quark-gluon content. 
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4. Behavior at large T 




Figure 6: High T behavior of the spatial masses of the ir and transverse p modes 
from the rank-1 separable model. 

The quark deconfinement point Td and the chiral restoration point T c are gener- 
ally expected be to identical or nearly so E3. (The present separable model produces 
Td = T c in rank-1 and Td = 0.9 T c in rank-2.) It might be expected therefore that 
above T c meson modes should dissolve in favor of a gas of essentially massless quarks. 
However for a significant temperature range above T c , the spatial tt and p modes 
studied here continue to be stable against qq dissociation and do not dissolve into 
a free quark gas. This situation is clearest for the simpler rank-1 separable inter- 
action. The results for M^{T) and Mj(T) obtained from the eigenmass condition 
A(— M 2 ) = 1 are displayed in Fig. |6| along with the quark dynamical mass function 
at p = 0. The masses of both spatial meson modes approach the asymptotic he= 
havior 2ttT from below. This asympt otic, b ehavior has been discussed previously 123 
and is observed in lattice simulations 

The manner in which the 2irT behavior at large T emerges from the present 
description is illustrated well by the rank-1 model. For T > T c , the dynamically 
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generated mass function of the quarks is essentially negligible, and the quark prop- 
agator is dominated by the Dirac vector amplitude uviq 2 ) ~ l/q n - For the spatial tt 
mode characterized by P = (0,P), the loop integral for the "polarization" function 
or BSE eigenvalue A,r(P 2 ), given by Eq. (||), yields at large T 



(2tt) 3 \k 2 T 2 + (q + P/2) 2 } [7r 2 T 2 + (g - P/2) 2 ] ' 

(35) 

where only the dominant zeroth fermion Matsubara mode has been retained. For 
T ^> Aq/tt, with Ao being the range of the interaction form factor /o, only small 
q is relevant and the position of the lowest singularity as a function of P 2 < ap- 
proaches 2ttT. The higher the temperature, the more A,r(P 2 ) is suppressed except 
near the singularity. The value A^— M 2 ) = 1 must be encountered before the di- 
vergence and thus M w (T) ~ 2irT — A w (T) where A n is a positive mass defect that 
will typically decrease with T. Thus the spatial meson mass or screening mass will 
approach the thermal mass of a pair of massless fermions from below. This limit 
has also been demonstrated from the pseudoscalar correlator within the Nambu- 
Jona-Lasinio model l!± 

In the general case, the detailed temperature behavior of the mass defect A 7r (T), 
or the nature of the approach to the 2-7rT limit, depends upon the asymptotic 
behavior of the quark amplitudes A(q 2 ) — 1 and B(q 2 ). In QCD, their leading 
asymptotic behavior is ~ 1/g 2 apart from slow logarithmic corrections. The form 
factors chosen in this initial work within a separable model induce an exponential 
fall-off. It is to be expected therefore that the present model estimate of A w (T) will 
decrease too rapidly with T. 

In Fig. the 7r and transverse p spatial masses at T > T c from both the rank- 
1 separable model and the ID model (see the Appendix for more details on the 
ID model) are presented in comparison with lattice QCD simulations of spatial 
screening masses Ej. The solid horizontal line marks 2tt while the lower horizontal 
dot-dashed line represents the lattice free limit corrected for the lattice time 
extent N t . From Fig. fj] it is evident that for the tt at T > T c , the spatial or screening 
mass defect A = 2irT — M(T) is decreasing more rapidly above T c in the rank-1 
separable model than is evident from the lattice simulations. This is consistent 
with the exponential behavior of the employed form factors. It is also consistent 
with the absence of quark vector self-energy amplitudes A(p, T) — 1 and C(p, T) — 1 
through which interactions can be quite persistent in the asymptotic region. This 
can be demonstrated within the chiral limit ID model where those amplitudes are 
strong and indeed have the power law fall-off. As seen in Fig. [?], the resulting mass 
defect for both 7r and p is in fact too strong when compared with the lattice QCD 
simulations. This persistent self-interaction well above T c , which slows the approach 
to free behavior such as Stefan-Boltzmann thermodynamics 113, may well be what 
is signaled by the lattice QCD data in Fig. |^. 



5. Discussion 
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Figure 7: Spatial masses from the rank-1 separable model (sep) and the in- 
frared-dominant model (ID) compared to spatial screening masses from lattice QCD 
simulations (latt) taken from E3. 



We have explored it and p spatial correlation modes at T > within the rainbow- 
ladder truncation of the quark Dyson-Schwinger equation and the qq Bethe-Salpeter 
equation in the Matsubara formalism. With parameters fitted to T = properties, 
the model possesses dynamical chiral symmetry breaking and quark confinement. A 
simple separable form of the effective interaction is employed and this facilitates the 
use of sufficient Matsubara modes to allow coverage of low temperatures as well as 
the transition region. Deconfinement and chiral restoration transition temperatures, 
Td and T c , are very similar and in the range 100-150 MeV. The T-evolution of 
the 7T and p qq states in the presence of the deconfinement and chiral restoration 
mechanisms is studied. The degree to which the model respects the axial vector 
Ward-Takahashi relation is evaluated in terms of the exact pion mass relation and 
the related GMOR relation. The p electromagnetic coupling constant g p and the 
strong coupling constant (7 p7r7r are also obtained as a function of T. Estimates are 
made for the T-dependence of the widths for p° — > e + e~ and p — > tttt. Finally, the 
high T behavior of the spatial masses is compared to that of spatial screening masses 
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from lattice QCD simulations c3. 

The masses M W (T), Mj(T) and Mp(T) are found to be almost T-independent 
below T c followed by a strong increase. This behavior is characteristic of lattice 
QCD simulations Ejtj and DSE studies Our estimate of the p — ► nir strong 

decay width shows a decrease with T such that, by T c , it has been reduced by 
50%. Our arguments suggest a phase-space blocking effect at about 25 MeV above 
T c . (A similar phase-space blocking has been argued before^jmly for the strong 
7T7T decay of the scalar-isoscalar partner of the pion near T c EJI23.) The tendency 
here is for the transverse p mode above T c to be left with a narrow total width 
typical of electromagnetic decay. One would expect the mass of the pseudoscalar 
K correlation to rise with T in a similar fashion to M T , while the masses of the 
vector <fi and K* modes should rise like the p. This suggests that the vector modes 
p, K* and tend to be trapped with their relatively long electroweak lifetimes 
and with significantly increased masses for a domain of high temperatures above 
the transition. Between 0.5-1.0 T c , the transition probability connecting vector 
correlations with pairs of pseudoscalars would be reduced. This suggests that within 
the gas of pions and other pseudoscalars that dominate the hot hadronic product 
from heavy-ion collisions, the role of vector meson correlations in producing the 
dilepton spectra could be significantly less than conventional expectations. The 
present findings follow from the response of the quark-gluon content of the mesons 
to the heat bath. A different phenomenon is the coupling of the meson modes to 
the many-hadron environment which introduces a collisional broadening effect □. 
An approach that incorporates both phenomena is clearly called for. 

Only the spatial or screening qq masses have been investigated within this model. 
The temporal masses (sometimes called dynamical or pole masses) provide a dif- 
ferent characterization of the correlations. Lattice simulations indicate that spatial 
masses become much larger than the temporal masses above T c EH. A Nambu- 
Jona-Lasinio model study c3 found that they are significantly different only in the 
range 150 MeV < T < 350 MeV. A possible explanation for this discrepancy lies in 
our finding that the mass defect A = 2irT — M(T) at high T is significantly influ- 
enced by residual non-perturbative interaction effects in the Dirac vector amplitude 
A(P 2 ) — 1 of the quark self-energy. Such a term is not present in the Nambu- Jona- 
Lasinio model which produces a momentum-independent quark constituent mass. 
The present model allows low temperature confinement and a momentum-dependent 
quark self-energy in a simple way but at a cost of an exponential asymptotic fall-off 
instead of a more realistic power law fall-off. Comparison with lattice QCD results 
illustrates the connection between asymptotic behavior of the interaction and the 
high T behavior of the mass defects A,r(T) and A p (T). 

The simplicity of a separable representation of the effective quark-quark inter- 
action of the type studied here might be of advantage in the consideration of pion 
loop effects and bulk thermodynamic properties of the hadron-quark matter ph ase 
transitions. Initial thermodynamic considerations have been reported recently cB. 
The only T-dependence of the effective quark-quark interaction implemented by 
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the present separable model is that generated by the Matsubara frequencies that 
enter through the momentum dependence. The strength and range parameters 
have been kept T-independent and no attempt has been made to introduce explicit 
T-dependent characteristics such as a Debye mass. Such considerations are more 
appropriately handled in approaches that have a better connection to a perturbative 
gluon propagator E3. 
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Appendix A 

The extension to T > of the ID model introduced at T — by Munczek and 
Nemirovsky El provides a semi-analytic perspective on the large T behavior of the 
meson masses. For the Feynman-like gauge used in the present study, the effective 
interaction of this model is specified by 

D(p-q)^{2n)^5\p- q ) , (A.l) 
16 

which is to. be used in the DSEs given in Eqs. (|) and (§) and also in the BSE given 
in Eq. (J^)Clln the chiral limit, closed form expressions exist for the resulting quark 
propagator amplitudes and for the ladder BSE eigenvalues A(P 2 ) in the pseudoscalar 
and vector channels of interest here. The chiral limit DSE solution is of the general 
form given in Eq. ( flcf ) with 




(A.2) 



TThe results are unchanged when Landau gauge is used if Eq. (^.1) is scaled up by 4/3. 
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The quark mass function is m 2 (p 2 ) — rj 2 /A — p 2 for p 2 < rj 2 /A 1 and m 2 (p 2 ) — oth- 
erwise. There is no mass-shell, i.e. p 2 + m 2 (p 2 ) ^ for any p 2 , and there is quark 
confinement. The above solution holds at T — and at T > with p 2 replaced by 
p 2 n = ijj 2 x +p 2 . Thus for T >T C = i]/(2ir) there is chiral restoration (gb = 0). Re- 
sults from this model at finite T, p have been obtained for quark thermodynamics 
outside the phase boundary of chiral restoration as governed by the quark propaga- 
tor E£l. The behavior of the masses of both the it and p modes for T <T C (andalso 
for chemical potential dependence for p < p c ) have been obtained previously E3. 



With Eq. ( |A.1|) , the BSE given in Eq. (|l|) becomes 



2 

X(P 2 )T(q; P) = -!L^S(q + )T(q; P)5(?_) 7/1 , (A.3) 

and solutions are possible if q is fixed by P. At T > 0, the solution that connects 
smoothly to the T = solution has the normally independent variable q n = [u) n , q) 
restricted to q n — ► (wo,0). To obtain spatial meson modes to compare with the 
separable model results, we again set P = (0, P). Within the temperature domain 
T <T C = T]/(2tt), the results for the pseudoscalar ancLvector spatial modes are par- 
ticularly simple and have been discussed previously tL3. One finds A^O) = 1 for the 
chiral limit it; thus — 0. For the vector meson one finds that Xj(— r] 2 /2) = 1; 
thus Mj = ri/^/2. These are also the correct T = results of the model. They hold 
over the finite temperature domain for which the quark mass function appropriate 
to the propagators occurring in the BS Eq. ( |A.3| ) is nonzero. For the equation 
appropriate to a meson of mass M, the relevant domain for the present model is 
rj 2 — As > where s — tt 2 T 2 — M 2 /4. For M„ = 0, this temperature domain corre- 
sponds to that for which the model DSE generates a dynamical quark mass function 
in accord with the Goldstonc theorem. For larger M, this temperature domain will 
be larger. The vector result M p = r\j\[2 holds for the larger temperature domain 
T < yfij2T c . We fix the single parameter r) = 1.107 GeV, so that M p = 0.783 GeV. 
This produces T c = 0.176 GeV. 

Beyond the temperature domain where the meson mass is constant, one finds 
A,r(P 2 ) = rf s_ <ta{s+) 2 , where s± — n 2 T 2 ± P/4, and the p eigenvalue is simply 
given by Xj(P 2 ) = ^Xt T (P 2 ). In general both functions A(— M 2 ) decrease with 
increasing T and increase with increasing M. Thus the spatial eigenmode condition 
X(-M 2 ) = 1 will tend to maintain Mj(T) > M„(T) while both masses rise with T, 
and the obtained behaviour of the masses Mj and M n is qualitatively the same as 
what we found in the separable model. The main difference is that the mass defect 
A(T) = 2ttT - M(T) is significantly larger in the ID model. The leading large T 
behavior 

Mj(T) 2 -> (2nT) 2 (l_ JL...) , (A.4) 

leads to the spatial mass defect of the transverse p mode having the asymptotic 
value A p (T) — > rj. Fig. fj] displays the behavior for both 7r and p. 
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